Non-Markovian Quantum Trajectories Versus Master Equations: 

Finite Temperature Heat Bath 

Ting YuJ] 

Rochester Theory Center for Optical Science and Engineering 

and 

Department of Physics and Astronomy, 
University of Rochester, 
Rochester, NY 14627, USA 

(Dated: 29 March 2004) 

Abstract 

The interrelationship between the non-Markovian stochastic Schrodinger equations and the cor- 
responding non-Markovian master equations is investigated in the finite temperature regimes. We 
show that the general finite temperature non-Markovian trajectories can be used to derive the 
corresponding non-Markovian master equations. A simple, yet important solvable example is the 
well-known damped harmonic oscillator model in which a harmonic oscillator is coupled to a finite 
temperature reservoir in the rotating wave approximation. The exact convolutionless master equa- 
tion for the damped harmonic oscillator is obtained by averaging the quantum trajectories relying 
upon no assumption of coupling strength or time scale. The master equation derived in this way 
automatically preserves the positivity, Hermiticity and unity. 

PACS numbers: 03.65.Yz, 42.50.Lc, 05.40.Jc 

Keywords: Open quantum systems; Non-Markovian quantum trajectories, Master equations 



* Also at Department of Physics, Queen Mary College, University of London, London, El 4NS, UK 
Electronic address: ting@pas.rochester.edu 



1 



I. INTRODUCTION 



Open systems are generic models for the study of quantum dynamics in the sense that a 
real quantum system is either difficult to be isolated from the influence of its environment 
or is deliberatelyput in touch with some purposely engineered devices in order to make a 
measurement In any case, the system and the environment, initially independent, will 

become entangled due to the interaction, then the system state will not remain in a pure 
state and the dynamics is described by a non-unitary process. In many physically relevant 
cases, a typical open quantum system normally involves a small system of interest coupled 
to a large system with large number of degrees of freedom, commonly known as heat bath, 
reservoir or generally environment. Traditionally, the dynamics of a system of interest is 
described by a Lindblad master equation which can often be derived if the standard Born- 



mm. 



Markov approximation is assumed Moreover, it turns out that such a Lindblad 

master equation also allows to be unravelled into various stochastic Schrodinger equations 
known as quantum trajectories Both the master equations of the 

Lindblad form and their stochastic unravelings have become an integral part of theories of 
open quantum systems. However, when the heat bath memory effects are relevant, such as 
in the cases of high-Q cavity, atom laser or structured environment, where the Born-Markov 
approximation ceases to be valid, hence the dynamics of the open quantum system must be 
described by a non-Markovian process Q| . It is lorig known that the derivation of a 
non-Markovian master equation is a formidable task jlfil . Iffij ] . 

Recent research on open quantum systems has suggested that, alternatively, the non- 
Markovian dynamics may be described by a diffusive stochastic Schrodinger equation known 



as non-Markovian quantum trajectories or quantum state diffusion equation 
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For an A^-dimensional Hilbert space, the stochastic 



Schrodinger equation evolves an A^-dimensional vector, so offers numerical advantages over 
a master equation which evolves an A^x density matrix. Very recently, we have shown that, 
beyond the numerical advantages and the conceptual merits, the non-Markovian stochastic 
Shrodinger equations may also provide a powerful tool for deriving the corresponding non- 
Markovian master equations. This idea has been explored in several distinct cases including 
a two-level atom interacting with a zero-temperature heat bath 2^ 28 1, and a Brownian 



particle coupled linearly to a finite temperature heat bath via the position variable commonly 
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known as quantum Brownian motion model 

0, y, Q- 

The purpose of this paper is to 
extend this research to more general finite temperature regimes where the Lindblad operator 
is not a Hermitian operator. In this paper, we take the damped harmonic oscillator as our 
primary example. This simple model is of great interest in quantum optics because it is an 
essential ingredient in the theoretical investigations of various quantum optical experiments. 
We show here in detail that the non-Markovian quantum trajectories allow the derivation of 
the exact convolutionless master equation irrespective of the coupling strength, the separated 
time scales, or the special distributions of the environmental frequencies. 

The organization of the paper is as follows. In Sec. m we introduce both zero and 
finite temperature stochastic Schrodinger equations. In Sec. 11111 we show how to derive 
an exact master equation from the finite temperature quantum trajectories. In Sec. IIV|, 
we establish the stochastic Schrodinger equation for the damped harmonic oscillator and 
show that the corresponding exact convolutionless master equation can be obtained by 
averaging the solutions to the non-Markovian stochastic Schrodinger equation without any 
approximations, in particular, without Markov approximation. We conclude this paper in 
Sec. El 



II. NON-MARKOVIAN STOCHASTIC SCHRODINGER EQUATION 
A. Zero-temperature (T = 0) 

Our model consists of a system of interest coupled linearly to a large number of har- 
monic oscillators with distributed eigenfrequencies uj\ and creation and annihilation opera- 
tors h\,h\. The quantum Hamiltonian for the system plus reservoir can be typically written 
as (we set h= 1) 

= H^y, + J2(9lL% + gxLb\) + ^Ah, (1) 

A A 

where g\ are the coupling constants and the system operator L coupled to the environment 
is often called the Lindblad operator. 

For the open quantum system described by IQ, the linear non-Markovian quantum state 
diffusion (QSD) equation has been derived from the formal solution of Schrodinger equation 
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for the total system in a special representation |18j : 

idt\^t) = Htot{mt), (2) 

where stands for the pure state vector for the total system and Htot(t) is the Hamiltonian 
Htat in the interaction representation with respect to the free bath Hamiltonian: 

X 

= ffsys + Yl {alL^^^xe-'-'' + gxLaie'^^'^ . (3) 

A 

In this subsection, let us assume that the initial pure state of the system and the environment 
is taken to be 

l^o) = l^o) ® |0i) ® IO2) • ■ ■ ® |0a) ® ■ ■ • , (4) 

with an arbitrary system state \iIjq) and the environment in the ground state |0). By us- 
ing a Bargmann coherent state basis for the environmental degrees of freedom: \z\) = 
exp{zxa\}\Ox) and the resolution of the identity 



e"l^^l>A)(^A|, (5) 



TC 



the total state |\E'() can be expressed as 



1**)= / — e-'^n^,(z*))®|z), (6) 



where l^;) = \zi) ® Z2) ® ■ ■ ■ ^ \zx) ■ ■ ■ , cPz = <Pzi<Pz2 ■ ■ ■ <Pzx ■ ■ ■ and \z\'^ = \zx\'^- Then 
the resultant pur e state for the system of interest {tptiz*)) was shown to satisfy the following 



equation 



;ant pur ( 
HQ: 



dtipt = -iHsys^pt + Lz*ilJt - dsa{t-s)-—^. (7) 

where a{t — s) = \gx\'^e^'^'^^^^^'^'^ is the bath correlation function, and z^ = —i g\z\e''^^^ 
is a colored, complex Gaussian noise with = }A[ztZs] = and ^^[2;^*^^] = a{t — s). 

Note here A^[-] is the statistical mean over the Gaussian process Zt- By construction, 

p, = MUt)m = I ^e-\^\'\i;,){U (8) 

The linear non-Markovian QSD equation ((7j) is valid for a zero-temperature heat bath, 
and also for a finite temperature heat bath with the condition L = . In general, if L 7^ L^, 
we shall see below that the finite temperature non-Markovian QSD equation takes a more 
complicated form. 



B. Finite temperature (T 7^ 0) 



The finite temperature stochastic Schrodinger equation or non-Markovian QSD equation 
can be obtained by mapping the system with the total Hamiltonian ^ and an initial thermal 
state to an extended system with a vacuum state such that the zero temperature stochastic 
Schrodinger equation for the extended system is equivalent to the situation where the heat 
bath is at finite temperature 



ISL Isi l- To be specific, for the total Hamiltonian we 
assume that the heat bath is in a thermal equilibrium state at temperature T, with density 
operator 

Pbath(O) = ^ , (9) 

where Z = Tr[e~^^'"'*] is the partition function and (3 = l/ksT. 

Now we introduce a fictitious heat bath with the bosonic operators c\, c\, which has no 
direct interaction with the system, ensuring that after tracing over the fictitious variables, the 
initial state of the original bath will be in the thermal state pbath(O). The total Hamiltonian 
with two independent heat baths is given by 

Htot = i^sys + Y,(9*xL^bx + gxLb{) + ^ ujxb{bx - ^a^ca- (10) 

A A A 

The boson Bogoliubov transformation formally couples the system of interest to two set 
of bosonic operators d\, d\ and e\,e\ : 



bx = Vnx + Idx + Vn^e{, (11) 
cx = Vnx + lex + v^4' (12) 

where nx is the mean thermal occupation number of quanta in mode Ux: 

exp{huJx/kBT) - 1 
The transformed Hamiltonian Ti'tot in terms of dx, ex is then given by 

n'tot = H,ys + Yl V^A + 1{9*xLUx + gxLd\) + ^ ujxd{dx 

X X 

+ \^{9xL'^A + 9\Lex) - ^ ^xe\ex. (14) 

A A 

Thus we have mapped the finite temperature problem into a zero temperature one with 
the initial vacuum state denoted by |0) = |0)d ® |0)e, satisfying dx\Q) = 0,eA|0) = 0. 



With the new Hamiltonian ()14p . the finite temperature problem has been reduced to a zero 
temperature one. Thus the resultant pure state for the system of interest ipt = \ipt{z*,w*)) 
satisfies the following stochastic Schrodinger equation with two independent noises z^,wl: 

Jo 

+L^w;^t-L [ dsa2{t-s)P^^, (15) 
Jo "^ly* 

-s) = J2^n, + l)\g,\'e-^-^^'-^\ (16) 

A 

a^it^s) = ^nA|^7A|V-^(*-^) (17) 



where 



A 

are the bath correlation functions and 



A 

< = 5^ ^?>Ie-*"^* (19) 

A 

are two independent, colored, complex Gaussian noises satisfying 

M[zt] = M[ztZs] = 0, M[z;zs] = - s), (20) 
M[wt] = M[wtWs] = 0, M[w;ws] = a2{t-s). (21) 

Note here A^[-] is the statistical mean over the Gaussian processes zl and wl- Again, 
by construction, pt = A4[\ipt){4't\]- In the zero temperature limit T ^ 0, we have ai{t — 
s) |(7aP exp(— iti^Al^ — s)) and a2(t — s) ^ 0, then equation (fT3j) reduces to the 

simple zero-temperature case ((Tj). However, in this paper, without explicit statement, we 
will always work in the finite temperature regimes. From equation (|T5|l . we see that the 
finite temperature heat bath has induced both the spontaneous transitions and stimulated 
transitions, moreover, it also gives rise to an absorptive process caused by taking thermal 
quanta from the heat bath. 

The stochastic Schrodinger equation (fT3j) can be greatly simplified by using the ansatz, 

^ = 0,{t,s,z*,w*)^t, (22) 

^ = 02{t,s,z*,w*)^t. (23) 



Then equation ()15|) takes a more compact form: 

+L^w;iJt-L02{t,z*,w*)iJt, 

where Oi {i = 1,2) denote 

Oi{t,z*,w*)= / ai{t — s)Oi{t, s, z* ,w*)ds, (i = l,2). 



(24) 



(25) 

Jo 

We can determine the operators Oi,2{t, s, z* ,w*) in (j^^ and (j^Hj) from the following 
'consistency conditions', 

5 



dt— = —dtiJt = —dt^Jt 

5zl 5zl ' 5wl 6w* ' 

together with the initial conditions: 

Oi{t = s, s, w*, z*) = L, 02{t = s, s, w*, z*) = 

From the consistency conditions, we may get the evolution equations for Oi 
Oi{t, s, z*,w*) and O2 = 02(t, s, z*, w*) 

dtOi = [-iH,y, + Lzt + L^wl - L^Oi - LO2, Oi] 



(26) 



(27) 



6zt 6z* 



(28) 



and 



dt02 = [-iH,y, + Lzl + L^wl - L^Oi - LO2, O2] 

-Ltfl-L^, (29) 

ow* ow* 

It should be remarked here that the 0-operators can be determined in many interesting 
situations and moreover, that the approximate 0-operators can always be obtained 



by invoking a perturbation technique 
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III. NON-MARKOVIAN MASTER EQUATION AT FINITE TEMPERATURE 

From the non-Markovian QSD equation ()24|) . one may take the statistical mean to derive 
the master equation, and gets 

dtpt = -t[H,ys,Pt] + [L, M{PtOl{t,z\w*)}]-[L\ M{Oi{t,z\w*)Pt}] 

+[Lt, M{PtOl{t,z*,w*)}] - [L, M{02it,z*,w*)Pt}]. (30) 
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This is the general master equation at finite temperature. Although this last result is still 
not a closed evolution equation for p^, it nevertheless shows how a convolutionless master 
equation may result from our knowledge of the operators Oj(t, s, z*, w*), {i = 1, 2). 

As we mentioned before, in the case of finite temperature with L ^ L\ the 0-operators 
will generally contain noises z^, w^. However, we begin with some special cases in which the 
operators Oi(t, s, z*, w*) turn out to be independent of the noises z^,w^. In such a case, we 
may write Oi{t, z*,w*) = Oi{t) and (jHUI) is indeed a convolutionless, closed master equation 



dtpt = -t[H,ys,Pt] + [L,PtOl{t)] + [Oi(t)pi,Lt] + [L\ptOlit)] + [02{t)pt,L] (31) 



for the non-Markovian open system. 

As a direct application of pip, let's discuss an important dephasing process: the system 
Hamiltonian commutes with the Lindblad operator [Hgys, L] = 0. More generally, we may 
consider a solvable model with [L,Hsys] = iK.I,L = = real constant, and / is the 

identity operator. For this case, it is easy to check that the solutions to ()28|1 and (jSHI) are 
given by 



By simply inserting and into (jHH), a closed convolutionless master equation is 
immediately obtained: 



-pt = -t[H,ys,Pt] + fmipuL] + r{t)[L,ptL]+g{t)[pt,L]+g*{t)[L,pt], (34) 



02(t, s, z* ,w*) = L — K(t — s). 



Oi(t, s, z* ,w*) = L — K(t — s), 



(32) 
(33) 



where the time-dependent coefficients are given by 




(35) 



Note here that the finite temperature correlation function a{t — s) is given by 



a{t — s) = ai{t — s) + a2it — s) = rj{t — s) + ivi^t — s). 



(36) 



where 




(37) 




(38) 



A 
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The master equation can be immediately applied to qubits decoherence and disen- 
tanglement [3^. To be specific, let's consider the dephasing process in a quantum register 
containing non-interacting qubits. It is known that one of the most important decoher- 
ence processes in the quantum registers is that qubits collectively interact with some 
environmental noise which only randomly interrupts the phase of each qubit. This problem 
can be effectively modeled by choosing i^sys = J2iLi ^i^t L = general 
Hamiltonian where al denotes the Pauli matrix for the ith qubit. The master equation 
for such a dephasing model is given by ()34j) with /t = 0. 

Of particular interest is the Markov approximation where ai{t—s) = 7i(5(t— s), a2(^— -s) = 
72(5(t — s), then Oi = ^L, O2 = -2■-^^ so the resulting master equation takes the well known 
Lindblad form 

+ ^ (2Ltp,L - Wpt - ptW) . (39) 

Before ending this section we have two remarks in order. First, let's note that the non- 
Markovian property of a convolutionless master equation ()31|) or more generally ()30p is 
characterized by the time-dependent coefficients. Second, if the 0-operators appearing in 
fl22|) and ()23p or in ()30|) do contain noises, as seen from the next section, things become much 
more complicated. There is no general guideline how to derive the closed convolutionless 
master equation from (j3(jp . However, it has been shown that the evolution equations for the 
0-operators Oj(t, s, z* , w*) with respect to s rather than t are essential for such a derivation. 
Here we believe that "The existence of a set of uncoupled evolution equations for the O- 
operators with respect to s" is a very powerful ansatz. Obviously, this ansatz remains to be 



tested in more examples 



29|. 



IV. THE DAMPED HARMONIC OSCILLATOR 



A. Stochastic Schrodinger equation for the damped harmonic oscillator 

In this section we show that the non-Markovian quantum QSD approach is versatile 
enough to handle the general finite temperature cases with L 7^ To simplify the cal- 
culations as much as possible, we consider here the simplest model of this kind but still 
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worth studying - the damped harmonic oscillator model. The model consists of a harmonic 
oscillator with frequency Q, coupled linearly to a large number of harmonic oscillators: 
i^sys = ^a^CL, L = a. We assume that the heat bath is at finite temperature T. This model 
has been discussed in various contexts (see, e.g., jsjO]). In what follows, we use this model 
to show explicitly how the finite temperature master equation can be derived directly from 
the non-Markovian quantum trajectory equation. For the damped harmonic oscillator at 
zero temperature, it has been shown that = and Oi = F{t)a, O2 = Q- Thus the 
zero-temperature non-Markovian master equation is immediately obtained from ()31|) . For 
the finite temperature heat bath, as to be seen below, the 0-operators Oi and O2 for the 
damped harmonic oscillator will contain noises and w^, hence things become much more 
involved. We emphasize that, for a quantum system with L ^ L\ the dependence of the 
0-operators on the noises is a generic feature for the finite temperature cases. First, note 
that the exact quantum state diffusion equation for the damped harmonic oscillator can be 
written as 

dti't = —i^a)aipt + o-zlipt — a)Oi{t, z* , w*)ipt 

+a^w;^lJt - a02(t, z*,w*)'ilJt, (40) 

where Oi are defined in 

Crucial to the solution of QSD equation ()4U|) is the explicit determination of the O- 
operators. It is easy to check that the following ansatz will give rise to the solutions of 
equation and 

Oi{t,s,z*,w*) = fi{t,s)a+ [ ds'ji{t,s,s')w;,, (41) 

02it,s,z*,w*) = f2it,s)a^+ [ ds'j2it,s,s')z:„ (42) 

^0 
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where fi{t,s) and ji{t, s, s') satisfy the coupled non/mear integro-differential equations: 

dthit, s) - znMt, s) - hit, s) f ds'aiit - s')h{t, s') 

Jo 

- hit, s) [ ds'a2{t - s)h{t, s') + f ds' a2{t' - s)j2(t, s' , s) = 0, (43) 
Jo Jo 

dth{t,s)+tnh{t,s) + h{t,s) [ ds' ai{t-s')h{t,s') 





+ h{t,s) [ ds' a2{t~s')h{t,s')+ [ ds"ai{t-s")Mt,s",s) = 0, (44) 
Jo Jo 

dtji{t, s, s') - hit, s) [ ds" ai{t - s")Mt, s", s') = 0, (45) 

dtj2it, s, s') + hit, s) [ ds" a2{t - s")32{t, s", s') = 0, (46) 
with the initial values: 

h{t = s,s) = l, h{t = s,s) = l, (47) 

ji{t = s,s,s') = 0, j2{t = s,s,s')=0, (48) 

ji{t,s,t) = -h{t,s), j2{t,s,t) = h{t,s). (49) 

Indeed we have seen from ()41|) and P2|l that, in the case of finite temperature, the O- 
operators generally involve the noises, even though at zero-temperature they do not. 

With the above evolution equations for hit, s) and ji(t, s, s'), we are able to simulate the 
system dynamics by solving the non-Markovian QSD equation pUll 3^ . But our aim in 
this paper is to take the quantum trajectories as theoretical tools rather than as numerical 
applications. Clearly, the solutions of the above nonlinear coupled equations are hard to 
deal with. In the next subsection, we will show remarkably that the evolution equations for 
hit, s) and hi(t, s, s') with respect to s can be turned into a set of first-order linear uncoupled 
equations. 

B. Decoupling evolution equations for O-operators 

The functions appearing in equations a s they stand, are very difficult to handle 

analytically. However, as shown in (also see \2^). great simplification may arise through 
investigating the dependence of the functions on s rather than t. It is proved here that the 
evolution equations for the O-operators with respect to s form an uncoupled set of linear 
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differential equations. The 
to the non-Markovian QSD 



cey observation comes from the Heisenberg operator approach 



2^ . For our purpose in this paper, the main task here is to 



find the evolution equations for 

A{s) = {z\{w\Utais)\0) = Oi{t,s,z\w*){z\{w\Ut\0), (50) 

and 

Bis) = {z\{w\Uta\s)\0) = O2it,s,z*,w*){z\{w\Ut\0), (51) 

where Ut is the unitary operator for the extended system (fT^ : \i'tot(t)) = l'^tot(O)). 
Clearly, the operators A{s) and B{s) depend on the time t, but for our purpose, we regard 
them, and therefore also the operators Oi(t, s, z*, w*), as functions of s. The time t appears 
as a parameter only. Obviously, the evolution equations for A{s) and B{s) are equivalent to 
that for Oi(t, s, z* ,w*) and 02(t, s, z* ,w*) with respect to s, respectively. First, note that 
the Heisenberg equation of motion for a(s) gives rise to 

tdMs) = ^ais) - J2 v/^^m gxe-'^'^'dxis) + J] gxe-'^^^'eiis). (52) 

A A 

Therefore, we get 

idsA{s)=nAis) + 5^v/im^?Ae-*"^^(^|(«;|WA(s)|0) 

A 

+ ^2-^^ 9xe~'^^^z\{w\Ute\is)\0). (53) 

A 

Simply integrating the Heisenberg equation of motion for the environmental annihilation 
operator d\{s), we get 



{z\{w\Utdx{s)\0) = -igWnx + 1 / ds' e'""'' {z\{w\Uta{s')\0), (54) 

Jo 

where we have used the fact that initially d\{0)\0) = dx\0) = 0. 

Next, we note that the operator e\{s) in (j^Hj) can be dealt with in a similar fashion. In 
order to find a closed equation, however, we have to carry out an expansion at the final time 
t rather than the initial value s = 0. This simplifies the calculations: {z\{w\Ute\(t)\0) = 
wl{z\{w\Ut\0) . Thus, we integrate the Heisenberg equation of motion for e\{s) given the 
final value at s = t to get 



e 



lis) = e{{t) - z^,gl I * ds' e-'-^''a{s'), (55) 
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hence, we get 

{z\{w\Ute\{sm=wl{z\{w\Ut\^)-i^>.g>. / ds' e-'^^' {z\{w\Uta{s'm. (56) 



Combining the results (jSU), and (j^Uj) . we obtain a hnear first-order differential 

equation for A{s)\ 

dA{s) + VlA{s) (57) 
-i / ds' ai{s - s')A{s') -i [ ds' a2{s' - s)A{s') = -iwl{z\{w\Ut\0), 







for s G [0, t] and a fixed final time t. Note that equation (j57j) has to be solved with the final 
value A{s = t) = a{z\{w\l(t\0). 

By noting that A{s) = Oi(t, s, z* , z*){z\{w\Ut\0) from (jnOj), so with the derived first 
order differential equation for A{s), we may immediately obtain the evolution equation with 
respect to s for the desired operator Oi{t, s, z* ,w*), we find 

idsOi{t, s, z* ,w*) = QOi{t, s, z* ,w*) — i / ds' ai{s — s')Oi{t, s' , z* ,w*) 

Jo 

—i / ds' a2{s' — s)Oi{t, s' , z* ,w*) + iwl, (58) 

J s 

with the final value 

Oi{t,s = t,z*,w*) = a. (59) 
Similarly, we get the evolution equation for 02(t, s, z* ,w*): 

ids02{t, s, z* ,w*) = —Q02{t,s,z*,w*) + i / ds' a2{s — s')02{t, s' , z* ,w*) 

Jo 

+i J ds' ai{s' — s)02{t, s' , z* ,w*) — iz*, (60) 

with the final value 

02it,s = t,z*,w*) = aK (61) 

Remarkably, we see that the coupled nonlinear equations (j^Hj) and ()29p have become a 
set of uncoupled linear equations (j^Hjl and (p?|) . We emphasize that equations and (jUnjl . 
with the final values ()59|) and respectively, are all that is required for the derivation of 
the master equation for the damped harmonic oscillator. 
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The equations (jSHj) and (lUUjl in terms of fi(t, s), f2(t, s) and hi{t, s, s'), h2(t, s, s') can be 
written as a set of uncoupled /m ear equations for s G [0,t], with t fixed, 

dJi{t,s)+tQf\{t,s)+ [ ds' a,{s-s')h{t,s')+ [ ds' a2{s' -s)j\{t,s') = 0, (62) 

Jo Js 

dsf2{t,s)-inf2{t,s)- ds' a2{s-s')f2{t,s')- ds' a,{s' - s) f2{t , s') = 0, (63) 

dsji{t,s,s')+tnj,{t,s,s')+ [ ds" ai{s- s")ji{t,s",s') (64) 

Jq 

+ j ds" a2{s' — s)ji{t, s", s') = S{s — s'), 

dsj2{t, s, s') - tnj2{t, s, s') - [ ds" a2{s - s")j2{t, s", s') (65) 

Jo 

- [ ds" ai{s" - s)j2{t, s", s') = -6{s - s'), 



with the final values 

/i(t,s = t) = l, f2{t,s = t) = l, (66) 
ji(t, t, s) = 0, j2(t, t, s) = 0, for all s. (67) 

C. Master equation for the damped harmonic oscillator 

Now we are going to derive the exact master equation for the damped harmonic oscillator. 
We will see that the evolution equations of Oj are crucial for deriving the corresponding non- 
Markovian master equation from the stochastic schrodinger equation (jlUj) . From (jHUjl . we 

get 

dtpt = -in[a^a,pt] + [a, M{PtOlit, z* ,w*)}] - [a^ , M{Oiit, z* ,w*)Pt}] 

+ [a\ M{PtOl{t,z*,w*)}] - [a, M{02{t,z*,w*)Pt}]. (68) 



For the equation (j68jl to be a closed, convolutionless master equation, we need the explicit 
expressions for the ensemble means: 

R,{t,s) = M[0,it,s,z\w*)Pt], R2{t,s)=M[02it,s,z\w*)Pt]. (69) 

For this purpose, we take the mean ■ ■ Pt] of equation (jSHj) . and then we get 

dsRi{t,s) +inRi{t,s) (70) 
+i [ ds' ai{s - s')Ri{t,s') -i [ ds' 02(3' - s)Ri{t, s') = -iM[wlPt]. 

Jo Js 
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Thus, we have to find an expression for Ai\ztPt]- By using Novikov's theorem 3^, it is 
possible to relate it to the 0-operator again j20|: 

M[wlPt] = [ ds'a;{s-s')M[PtOl{t,s',z\w*)]= [ ds'a;{s ~ s')Rl{t,s'). (71) 
Jo Jo 

Thus, we find the equations for the operator Ri(t, s), as function of s, to be 

dsRi{t, s) + inRi{t, s) (72) 
+ [ ds' ai{s - s')Ri{t,s') + [ ds' a2{s' - s)Ri{t,s') = [ ds'a^it - s')Rl{t, s'). 

Jo J s Jo 

and similarly, for the operator R2{t, s), we have 

dsR2{t,s) - inR2{t,s)- [ ds' 02(3- s')R2it,s') (73) 

Jo 

- ds' ai{s' - s)R2{t, s') = - ds'a\{s - s')R\{t, s'). 

Note that the equations (j72|) and (j73|) have to be solved with "final values" 

Ri{t, s)\s=t = apt, R2{t, s)\s=t = a)pu (74) 

as one can easily see from (fS^ and (jHH)- Note that equations (f7^ and (fTSj) are a set of 
coupled hnear equations that are expected from (PT| and (jl^ . 

By observing PT|) and (ji^ . we find that the solutions of the crucial equations (f?^ and 
(j73|) with the final values (j74|) must take the following form: 

Ri{t,s) = F{t,s)apt + G{t,s)pta, (75) 
R2it,s) = H{t,s)a^pt + Iit,s)pta\ (76) 

where F{t, s),G{t, s), H{t, s), and I{t,s) are complex functions with the appropriate final 
values at s = t: 

F{t,s = t) = 1, G{t,s = t) = 0, (77) 
H{t,s = t) = 1, l{t,s = t) = 0. (78) 

Once we have Ri{t,s) = M[Oi{t, s, z* ,w*)Pt] {i = 1,2)) in the form of ^ and dZBD, it 
is straightforward to write down the master equation. Simply inserting result (j75p and (|76|) 
into the general form of the non-Markovian QSD- master equation ()68|) . we get 

dtpt = -ifi[a"''a,pj] + a(t)[a, pa^] + b{t)[a, a^] + c{t)[a\ ap] + d{t)[a\ pa], (79) 
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with 



a(t) 



bit) 



ds [alit 
a*{t), 




(80) 



(81) 
(82) 
(83) 



c(t) 



d{t) 



b%t). 



The equation (f7^ is vahd for arbitrary temperature, so it is expected to provide a good 
description for the low temperature regimes where the Markov approximation is doomed 
to fail. Clearly, the master equation ()79|) derived in this way automatically preserves the 
positivity, Hermiticity and trace. Moreover, it can be easily shown that the above master 
equation is of the Lindblad form, but the coefficients are time-dependent, so it depicts the 
non-Markovian dynamics. 

D. Determination of the coefficients of master equation 

In what follows, we show that the coefficients a{t) , b{t) , c{t) and d{t) can be expressed 
in terms of solutions of some basic equations. Suppose that u(t, s) is the solution to the 
following homogeneous equation: 



Then from (f?^ and (fTSj) it can be shown that F{t, s), G(t, s), H{t, s) and I(t, s) satisfy the 
following inhomogeneous equation, 




(84) 



with the final value u{t,s 



t) = 1. Here the kernel function f3{s, s') is defined as 




(85) 




(89) 



(87) 



(86) 



(88) 
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where the functions X{t, s) and Y{t, s) are given by 

X{t,s) = [ ds'a2{s' - s)u*{t,s'), (90) 
Jo 

Y{t,s) = I ds'aiis' - s)u{t,s'). (91) 
Jo 

Hence it is easy to see that the coefficients of the master equation (fTUj) can be expressed in 
terms of the solutions of those basic equations (j8l|l and (j86|l -f 



V. CONCLUSIONS 



The non-Markovian quantum trajectories offer a brand new method for exploring a quan- 
tum system coupled to a non-Markovian environment. Such a situation appears in various 
physical problems, such as materials with photonic band gaps or output of atom lasers. It 
is known that master equations and quantum trajectories in Markov regimes are of funda- 
mental importance for the description of open system dynamics. Moreover, they are often 
complementary to each other to provide a full picture of the underlying physics. In this pa- 
per we show that in the case of general finite temperature heat bath this fruitful interrelation 
between the master equations and the quantum trajectories can also be established in the 
non-Markovian regimes. In particular, we show that by averaging the stochastic Schrodinger 
equation the exact non-Markovian master equation of the damped harmonic oscillator at 
finite temperature can be obtained. The master equation derived in this way takes the 
convolutionless form, so the non-Markovian property is encoded in the time-dependent co- 
efficients. 

In the Markov regimes, the derivation of a master equation poses no special difficulties. 
However, it is a rather difficult task to derive a convolutionless exact master equation in 
the non-Markovian regimes. As shown in this paper, non-Markovian quantum trajectories 
provide a very useful tool in handling exact or approximate non-Markovian master equations. 
We believe that the techniques employed in this paper may find broader applications to other 
quantum open systems, leading to the establishment of a general master equation for the 
system interacting with a bosonic heat bath. This will be the topic of future investigations. 
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